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Abstract 

We study effects of direct interatomic interaction on cooperative processes in atom-photon dy- 

. namics. Using a model of two-level atoms with Ising-type interaction as an example, it is demon- 

O _ 

CN I strated that interparticle interaction combined with atom-field coupling can introduce additional 

interatomic correlations acting as a phase synchronizing factor. For the case of weakly interacting 
atoms with J < fkVQ, where J is the interparticle coupling constant and ujq is the atomic frequency, 
. . dynamical regimes of cooperative relaxation of atoms are analyzed in the Born-Markov approx- 

I imation both numerically and using the mean field approximation. We show that interparticle 

' correlations induced by the direct interaction result in inhibition of incoherent spontaneous decay 

o : 

\^ . leading to the regime of collective pulse relaxation which differs from superradiance in nature. For 



o 



superradiant transition, the synchronizing effect of interatomic interaction is found to manifest 
itself in enhancement of superradiance. When the interaction is strong and J > Hujq, one-particle 
■ one-photon transitions are excluded and transition to the regime of multiphoton relaxation occurs. 

^3 : 

Q^, Using a simple model of two atoms in a high-Q single mode cavity we show that such transition 

^ is accompanied by Rabi oscillations involving many-atom multiphoton states. Dephasing effect of 

^ • dipole-dipole interaction and solitonic mechanism of relaxation are discussed. 

PACS numbers: 42.50Fx; 78.67.-n 
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I. INTRODUCTION 



The phenomenon of superradiance has a long history dating back more than 60 years to 
the seminal paper by Dicke [l| where the effect was predicted theoretically. Over the past 
few decades the superradiance has been the subject of intense theoretical and exgerjmental 
studies in a large variety of systems. These include molecular aggre gate s 



atoms ja, |7| and Bose- Einstein condensates |8 



nanoclusters 



li 



15 



3 



heterostructures 



J9, 

18, 



19 



IC I , atomic nuclei 



igate s m , Is 



cperin 

M 



cold 
magnetic 



and many others. 



The key process underlying the mechanism of superradiance is phase synchronization of 
initially independent atoms caused by the coupling with a common environment represented 
by the electromagnetic field. In order for such process to occur the phase decoherence time 



of atoms should be longer than the photon travel time in the sample 

For samples which size is smaller than the wavelength of radiation, this condition requires 
the density of atoms to be sufficiently high. The system-environment (atom-field) coupling 
also manifests itself as an additional indirect interaction (a sort of the transverse dipole- 
dipole interaction) which may suppress superradiant transitions depending on the spatial 
distribution of atoms or the sample geometry 0, |21|. 

From the other hand, when the density of atoms (or, more generally, emitters) is high the 
direct interparticle interaction starts to play an increasingly important part in determining 
cooperative behavior of the particles. In particular, this interaction strongly affects the 
properties of low-dimensional systems. The Mott-insulator quantum phase transition in 



optically trapped atomic systems 



22 



y, 24, 



and in solid structures 



generation of many-particle entangled states or many-particle coherent dynamics, as it i s 
in the case of effectively interacting atoms inside a high quality dissipative cavit 




Bose-Einstein condensate 23l . l32| and in molecular clusters with strong magnetic 

are examples. 

The direct interparticle interaction introduces additional correlations between emitters. 
These correlations considerably influence the cooperative optical properties of atoms. 

Firstly, the interaction directly affects the superradiance leading to a number of pecu- 
liarities such as changing the order of superradiant phase transitions 



37 



38l |. Recently, the 



possibility of superradiant relaxation in strongly correlated systems was studied theoreti- 
cally The experimental results for magnetic molecules of Mni2-ac type were also 
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reported in Ref. Q|. Earlier, systems of ferroelectric type with strong interparticle inter- 
action were regarded as promising candidates for an active medium of the heat pumping 
laser bs]. 

Secondly, the direct interparticle interaction can play the role of a phase synchronizing 
factor that may lead to the cooperative behavior which, though shares many common prop- 
erties with the superradiance effect, essentially differs from superradiance in nature. The 
classical example furnishes the spectrum of P luminescence band in CdS and ZnO where the 
emission intensity is proportional to the second power of the free exciton number (pumping 
inten^ty). In tMs case the effect i. caused by exciton-excton scattenng [1 Recently 
such effects were observed in the microcrystalline phase of CsPbCla thin films |41| . 

For interacting atoms, the interaction can drastically change the regime of atom-photon 
dynamics by inducing (otherwise, excluded) multiphoton transitions 4^. It was shown in 
Ref. j43| that interatomic interaction can give rise to non-zero multiphoton emission observed 
with single-molecule spectroscopy technique as a two-photon cooperative effect for strongly 
dipole-dipole coupled molecules. Theoretically, this phenomenon was predicted as a large 
two- atom two-photon resonant effect for two atoms inside a high-quality cavity 44] . 

So, different regimes of radiative decay in correlated atomic systems are mainly governed 
by the interatomic interaction. By controlling the interaction radiation properties of such 
systems can be widely varied ranging from superradiant transitions to the generation of the 
Fock state of light. In particular, such control is feasible for the atoms in optical lattices 
(see, e.g., HQ)- 

In this work cooperative radiation of interacting atoms coupled to electromagnetic bath 
will be of our primary interest. We are aimed to study different relaxation regimes deter- 
mined by the intensity of interatomic coupling. 

The paper is organized as follows. 

In Sec. |n]we formulate the model of two-level atoms with Ising-type interaction and 
qualitatively discuss various regimes of relaxation by considering realignment of the atomic 
energy spectrum at different values of the interatomic coupling constant J. There are two 
limiting cases of weak and strong interaction with J < TujJq and J > hwQ, respectively {ujq is 
the atomic frequency). 

We find that, for weakly interacting atoms, Ising interaction would affect dynamical be- 
havior of the system leading to the transition to collective pulse relaxation and enhancement 
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of super radiance. For strong interaction, the regime of multiphoton relaxation is predicted 
to occur. It is also shown that dependence of the radiation intensity peak on the number of 
particles can be anomalous at long-range interatomic interaction. 

Derivation of the master equation for weakly interacting atoms is presented in Sec. IIIII 
We show that, for certain atomic configurations, dephasing effects of induced dipole-dipole 
interaction can be suppressed and dynamics of the atomic system can be described by the 
simplified master equation. 

In Sec. II VI the effects for weakly interacting atoms briefly discussed in Sec. |n]are investi- 
gated in detail. By applying the mean field approximation we obtain the results that agree 
very well with those calculated by solving the equations for atomic variables numerically. 

The regime of multiphoton relaxation that takes place at strong interatomic interaction 
due to inhibition of one-particle one-photon transitions is described in Sec. Finally, in 
Sec. IVIl we draw together the results, discuss solitonic mechanism of relaxation and make 
some concluding remarks. 

Details on some technical results are relegated to Appendixes El- O 



II. ATOMIC ENERGY SPECTRUM AND REGIMES OF RELAXATION 



In order to illustrate a possibility of different relaxational regimes, caused by direct inter- 
atomic interaction, we consider a simplest model of a chain of two-level atoms with nearest 
neighbor Ising-type interaction. Typically, the models of this type are used to describe 
molecular systems with Coulomb and magnetic interactions such as ferroelectric 3| and 
magnetic j^, 48, 4^ clusters, interacting electrons in the tightly-binding approximation j^l, 
interacting atoms in optical lattice j^l, where the extended Hubbard model is usually in- 
troduced. 

The interaction takes into account repulsion of the neighboring atoms so that the system 
may reveal the correlated many-particle behavior. When the energies of the neighbors are 
different tunneling transitions between neighboring sites are inelastic. 

In addition to the Ising-type direct interaction, the atoms interact with a common electro- 
magnetic field (the atom-field interaction). These two interactions are combined to maintain 
correlations in the atomic subsystem thus crucially affecting the regimes of cooperative op- 
tical transitions. 
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A. The model 



The Hamiltonian of the model atomic chain with the short-range Ising-type interaction 
can be written in the following form 

H = Ha + Hf + Hi, (1) 

where 

N N 

HA = hwoJ2st-jJ2stS!^, (2) 

i=l i=l 

is the Hamiltonian of atomic subsystem, 

Hp = ^ ^fcflk^aks (3) 

ks 

is the Hamiltonian of electromagnetic field and 

N 

Hj = -ihJ2Yl bks, {St + S;) ak, - H.c] (4) 

i=l ks 

is the Hamiltonian of the atom-field interaction and H.c. stands for Hermitian conjugation; 
ujo is the frequency of atomic transition, J > is the coupling constant of the nearest neigh- 
bor interaction, is the number of atoms, S"? is the 2;-component of the pseudo-spin operator 
determining the population of the ith atom, (•S'j^) is the raising (lowering) operator of 
the pseudo-spin, a^s {(^^.s) ^^^^ photon annihilation (creation) operator characterized by 
the wavenumber vector k, the frequency Uk and the polarization vector e^s, 

is the coupling constant of the dipole interaction between the ith atom and transverse 
electromagnetic field, dj is the effective dipole moment of ith two-level atom with the vector 
of spatial coordinates rj, Eq is the vacuum dielectric constant and v is the volume of field 
quantization. 

The Hamiltonian is distinguished from standard models of two-level atom interacting 
with the electromagnetic field by the presence of the interatomic interaction in Ha [see 
Eq. ©]. 

In what follows the dynamics of the collective relaxation of initially inverted atomic 
subsystem will be of our primary concern. More specifically, we shall study various dynamical 
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scenarios that occur at different values of the couphng constant J. In the remaining part 
of this section we first discuss energy spectrum of the atomic subsystem and quahtatively 
describe these dynamical regimes. 



B. Weak interatomic interaction, J < fkjQ 

The transition energy of an atom depends on the states of its neighbors as a result of 
interatomic interaction. This effect can be described as renormalization of the transition 
frequency of the ith atom ujq Coi which can also be readily seen from the mean field 
expression for the Hamiltonian Ha- 



Sufficiently weak interaction, J <^ hujQ, results in inhomogeneous broadening of the ra- 
diation spectrum. When the coupling constant of interatomic interaction J increases and 
J < hujo, the spectrum of the atomic subsystem undergoes more pronounced rearrangement. 
As is shown in Fig. the spectrum is no longer equidistant, but its structure is identical to 
that for non-interacting atoms with J = 0. In particular, the latter implies that the state 
with all the atoms excited | T) = I T)T) •••)?) (spin up represents the excited state of the 
atom) gives the level of highest energy and the level of lowest energy corresponds to the case 
where atoms are all in the ground states | |) = | |,|,...,|). 

Similar to the case of non-interacting atoms, the dynamics of radiative decay is domi- 
nated by one-particle one-photon transitions and can be described using the Born-Markov 
approximation. 

A cooperative relaxation of the system from an excited state can be realized by different 
schemes of transitions. It is a possibility to have a generation in the two mode regime with 
the frequencies uj = ujq± J/h (see Fig. ^) arising due to renormalization of the atomic 
frequency uoi that governs the intensity of radiation. 



It is well known |21l] that the intensity of atom radiation / is proportional to the fourth 
power of the transition frequency / oc cugldp, where d is the transverse dipole moment of 
the atomic transition. For weakly interacting atoms, the intensity / is defined by the 
renormalized frequency and can be represented by a sum of the coherent and the incoherent 



J 
2h 




parts, / = Jcoh + /incoh, where hoh oc Xli^j ^^i'ldjd*] oc uo'^N^ and Jincoh oc ufdj oc uj' 
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FIG. 1: Energy spectrum of Ha for the chain of six atoms. Two cases are shown: (a) weak 
interatomic interaction with J/fkVQ = 0.1 and (b) strong interatomic interaction with J/Huq = 10. 
The state | T) = I T)T)"-)T) (I i) = I indicates that the atoms are all in the excited 

(ground) state. 

Dependence of the renormalized frequency on the states of neighboring atoms introduces 
additional correlations into the incoherent part of radiation. Such correlations may result 
in a cooperative radiation of fundamentally different origin than superradiance described by 
the coherent part Jcoh- 

The effect of superradiance, for its part, can be enhanced by the Ising-type interaction 
which may act as a phase synchronizing factor. It is also worth noting that, for a long-range 
interparticle interaction, the renormalized frequency Ui depends on the number of atoms, 
oji ~ uq (^1 — I Xlj^i ■^('S'j ) j 5 leading to anomaly in the iV-dependence of the intensity 
peak of collective radiation. 

Figure^ shows that, when the atomic subsystem is not completely inverted, an excited 
state may relax through a series of transitions which frequencies are identically equal to the 
resonant frequency Uq. This regime of relaxation resembles the solitonic mechanism where 
the process of emission is connected with inelastic motion of defects (Bloch walls) induced 
by a short-range interaction. We shall discuss this point at greater length in Sec. IVII 
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C. Strong interatomic interaction, J > fkvo 



In the case of strong interaction with J > hio^, arrangement of the energy level signif- 
icantly differs from that of non- interacting atoms. Now the state | f) does not give the 
level of highest energy and < E^^i, where Ej is the energy of the state with j excited 
atoms. Such realignment indicates that many-atom multiphoton transitions will dominate 
the process of relaxation. 

Figure [TJd illustrates that, in the limiting case with J ^ hujQ, the states | t) and | |) give 
two lowest lying energy levels. They are separated from the other excited states by the gap 
of the width ~ J. Relaxation of fully inverted state (transition from | |) to | |)) can be 
achieved through a A^-photon process. 

At this stage it should be stressed that an accurate description of the interaction between 
the system and the field in the presence of strong interatomic correlations requires careful 
consideration of interplay between many-atom multiphoton dipole and multipole transitions. 
In a sense, our system can be regarded as an intermediate case between a collection of two- 
level atoms independently interacting with the field and a single many-level system. 

A good example provides the spin-phonon interaction in magnetic molecules Mni20i2 
of the multipole form Hj oc {S~^)'^ak + (S'~)^a^ which symmetry is consistent with the 
anisotropy energy [s^ . Similarly, the spin superradiance in atomic nuclei can be dominated 
by multipole transitions Il3l l. 

n 

For the multipole superradiance , the exponent of the A^-dependence of the radiation 
intensity peak (A^ is the number of emitters) can be greater than two, / oc A^" with a > 
2. This is the case for the cooperative spin-phonon relaxation in Mni20i2 where a rough 
estimate gives I oc {{3+^(3-^) oc A^^. 

The standard approach to the dissipative dynamics of atoms relies on the Lindblad equa- 
tion for the reduced density matrix of the atomic subsystem derived by eliminating boson 
variables from the master equation for the density operator of atom-field system. This 
master equation is obtained using the Born-Markov approximation 52] that neglects mul- 
tiphoton transitions and thus is applicable only for weakly interacting atoms at J < Hooq. 
In the subsequent section we discuss the derivation procedure for the Lindblad equation for 
the density matrix of weakly interacting atoms. 
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III. MASTER EQUATION FOR WEAKLY INTERACTING ATOMS 



In this section we derive and discuss the master equation for the reduced atomic density 
operator p for the case of weak interaction at J < hwQ. We show that the direct interaction 
renormahzes both the induced dipole-dipole interaction and the damping operator in the 
Lindblad equation. In addition, we find that, under certain conditions, destructive effect of 
the dipole-dipole interaction on supperradiance is suppressed. 



A. Derivation of master equation 

We begin by applying the standard procedure of elimination of field variables in the 
Born-Markov approximation to the equation for the density matrix of the combined 
field- atom system, paf-, written in the representation of interaction 

ihj^PAF{t) = [V{t)rPAF{t)], (6) 

where 

pj,p{t) = e'iH^+^^)t/Hp^^^-^iH^+H,)t/H^ 

V{t) = e^iHA+Hp)t/hj^^^-i{HA+Hp)t/h ^ 

N 

i=l ks 

f, = 1 - /3 {St^, + SU) , (3 = J/{hwo) (8) 



and the expression for the operator ^ is derived using the cyclic boundary conditions. 

As it was pointed out in Sec. m for the ith atom, the energy of transition depends on the 
state of the neighboring atoms. Alignment of the energy spectrum of the atomic subsystem 
and, as a result, the character of atom- field interaction are determined by the value of the 
ratio P = J/hjjjQ. When J > hujo and /3 > 1, ordering of the energy levels changes and 
certain transitions appear to be excluded. 

In Eq. ((Zj) the frequency of the dipole transition is effectively described by the operator (jHl) 
uJoTi entering the Hamiltonian of interaction V. As opposed to the case with f3 < 1, the 
eigenvalues of the operator f j can be negative at /5 > 1 depending on P and i. 

In this case, for example, the energy of the "monodomain" state | j, j, |) is below 
the energy -Eat^i of the state where only one atom is in the ground state | |, |, f, .., f). 



The result is that relaxation of the completely inverted state | |, |, |) cannot occur as a 
one-photon process. 

It should be noted that, for (3 > 1, the change of sign in the eigenvalues of the operator Fj 
introduces new resonant terms of the form S^a^^. However, as it was discussed in Sec. Ill C\ 
this case is additionally complicated by multipole transitions that need to be taken into 
consideration. 

When J < fvjjQ and /3 < 1, the eigenvalues of the operator f j are all positive. So the Ising- 
type interaction just renormalizes the transition frequencies inducing additional correlations 
in the atomic subsystem. 

Dynamics of the atomic relaxation is now governed by one-photon processes and, in the 
remaining part of this section, we restrict ourselves to this case of weakly interacting atoms 
with (3 <l. 

We can now use the Born approximation to derive the equation for the reduced density 
matrix of the atomic subsystem [3], P = ^^fIpaf}- For simplicity, we consider the zero- 
temperature case where only vacuum fluctuations of the electromagnetic field are taken into 
account. So, the density operator paf can be written in the following form 

PAF{t) = \0){0\^p{t). (9) 

We additionally assume that the wavelength of radiation is much longer than the size of 
the system and use the Markov approximation [2^. Omitting the technical details that can 
be found in Appendix |X1 the resulting equation is given by 

7 poo f •'T^ 

"J TT Jo \ 00 — ooni. i 



+ [R,{t)rp{t)S-{t)^^] + H.c. \, (10) 



where 



and 



p{t) = eiiHA+Hp)t/hp^-i{HA+Hp)t/h^ ^-^-^s^ 

R,{t) = snt) + s;{t), 



= ^ {[1 - (d ■ + [1 - 3(d . r,)1 " ^ ) ^ (12) 
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k = uj/c, Tij = \rij\ = |rj — rj\ is the separation distance between ith and jthe atoms, we 
assume dj = d^ = d, fij = rjj/|rjj|, d = d/|d| and 

For brevity, in Eq. (fTUI) we used operator functions defined on eigenvalues of Fj, so that if 
ti\ip) = Ti\ip) then, e.g., 6{uj — c^of = \ip)S{uj — cuoFj). It should be stressed once again 
that Eq. (|Tn|l is valid only for the case of weakly interacting atoms, when J < hwo and the 
Born-Markov approximation is applicable. 

Principal values of the integrals that enter the right hand side of the master equation ()10|) 
define the dipole-dipole interatomic interaction (see Appendix El and the Bethe-part of the 
Lamb shift, which depend on the direct interparticle interaction. Equation ()10|) can be 
further simplified under the condition that the characteristic wavelength of radiation A is 
much longer than the size of the sample. 

It is well known that in this case the dipole-dipole interaction is of quasi-static char- 
acter and is independent of the frequency of atomic transition. Since the direct interaction 
in Eq. (|7j) just renormalizes the atomic frequencies ujq ^ Ui = ujQTi, the asymptotic form of 
the dipole-dipole interaction in the limit r^/X <^ 1 must be identical to that for the system 
of non-interacting atoms with J = 0. At the same time the relaxation terms proportional 
to ujf essentially depend on the interatomic interaction. 

Mathematically, Eq. ()A6j) in the limit rij/X <^ 1 gives the relation 

With Fij{uj) 1 and the Lamb shift disregarded Eq. (fTIH) gives the master equation in the 
simplified form 

j^P{t) = 70 Yy^AtlpmStit)] + ^ Y.^[R,{t)rp{t)Ut)] + (15) 

where 70 = t^tjuJq is one-half the spontaneous decay rate for an isolated atom and 

3 l-3(d.ri,)' 
= -"'2 ■ (r,/c)3 



is the constant of the quasi-static dipole-dipole interaction 

The first term on the right hand side of Eq. (fT3j) describes renormalization of the damping 
term by the interatomic interaction. After neglecting the rapidly oscillating terms (the 
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rotating wave approximation) in Eq. (fT3j) and going back to the Schrodinger representation, 
we derive the master equation for the case of weakly interacting atoms, /? < 1, in the final 
form: 

+ 70$^ ([^-,p(t)f?5+] + [S-t^p{t),S;]) . (17) 

id 

Eq. (|T7|l describes dynamics of the Ising-like interacting atoms coupled to electromagnetic 
bath. The key feature of this equation is that only the damping terms are renormalized by 
the direct interatomic interaction at rj,- ^ A. 



B. Effects of dipole-dipole interaction 

It is known (see, e.g., reviews 21 1) that the dipole-dipole interaction may destroy 
the cooperative radiation due to the frequency chirping effect. This effect is caused by the 
dynamical shift of energy levels that results in a rapid loss of phase synchronization with 



the characteristic time of phase decoherence much shorter than the rate of relaxation 

Influence of the dipole-dipole interaction on the dynamics of the atomic subsystem is 
strongly affected by the spatial configuration of atoms and the angular distribution of dipole 
moments. It turns out that, for certain configurations, dephasing induced by the dipole- 



dipole interaction is inhibited and the cooperative radiation takes place |2G , 

In particular, the latter is the case when atoms are all in identical environment and the 
system is invariant under permutations of atoms, Qij = Qi>j> = Q. These conditions are met 
for a two atom system or an atomic ring [2^. 

For disordered atomic configuration (atomic cloud), the superradiance effect depends on 
the sample shape. It can occur for pencil-shape patterns as a result of suppression of the 
destructive effect of the dipole-dipole interaction 21 1. 

In our case, when the interatomic spacings is much smaller than the wavelength (r^j ^ A), 
the form of the dipole-dipole interaction ()l(j|l is the same for both interacting and non- 
interacting atoms. Symmetry of the renormalized damping terms is also identical to the 
symmetry of the system of non- interacting atoms with J = 0. So, atomic configurations 
suppressing the dephasing effect of the dipole-dipole interaction are the same for both cases: 
J = and J ^ 0. 
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For the simplest case of two atomic configuration, this point can be illustrated explic- 
itly. Similar to the case of two non-interacting atoms considered in Ref. js^, the master 
equation ()17|1 for two atoms with J 7^ is exactly solvable. The exact solution described in 
Appendix El shows that relaxation is independent of the dipole-dipole interaction. 

The analytical results presented in Appendix El can also be used to demonstrate the 
effect of additional correlations induced by the interaction. For this purpose, we consider 
the expression for the relaxation rate of initially excited atoms 



l{t)^-jTr{p{Sl + S^,)} 



Ame 



-Amt 



1 + 



n 



n — m 



(1 



-i(n-m)t 



(18) 



where n = 70(1 + (3/2)^ and m = 70(1 — /?/2)^. 

In Fig. 121 the relaxation rate (fT^ is plotted as a function of time at various values of 
/3. It is seen that the interaction noticeably influences the character of relaxation even at 
small values of j3. By contrast to the case of non-interacting atoms with /3 = 0, the system 
features collective decay with a pronounced emitted pulse arising as a result of additional 
correlations induced by Ising-type interaction. 
4 




FIG. 2: Relaxation rate ^(t) ()18() for two interacting two-level atoms at various values of the ratio 
f3 = J/hwo: (a) /? = 0, (b) /? = 0.1, and (c) [5 = 0.2. 



Now it is our primary task to examine how the cooperative relaxation is influenced by the 
interatomic interaction. So, in what follows we shall restrict our study to the atomic con- 
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figurations where the frequency chirping effect can be disregarded. For such configurations, 
the master equation (fT7j) can be considerably simplified by neglecting the terms describing 
the dipole-dipole interaction. 

IV. COOPERATIVE RADIATION OF WEAKLY INTERACTING ATOMS 

From the master equation ()17|) without the dipole-dipole interaction, Qij = 0, it 
is straightforward to deduce a system of equations for averages of the atomic variables 
{S^} = TrA{pS^} and {S^) = TrA{pS^}. The result is given by 

j^iSn) = -^oJ2(^ISlS~ + S^Srrl), (19a) 

i 

= zo.o(^+f,„) +270 J](ff^+^^), (19b) 

i 

j^{S-) = -zuo{S-h) + 2^oYliS:.S-rl), (19c) 

i 

where 70 = nrjUQ is one-half the spontaneous decay rate for an isolated atom. 
The relaxation rate of the atomic subsystem 

= - E ^ = ^0 J2(^!S^S- + S^Srr^) (20) 

n n,i 

is of primary importance in our subsequent analysis. The parameter 7(t) is defined in the 
right hand side of Eq. ()19ap . 

By contrast to the case of non-interacting atoms, there is no an intimate connection 
between the relaxation rate 'y{t) and the total intensity of radiation I(t) which is proportional 
to J^ini^n^i^'^i + ^i"^i^n)- The reasou is that the interatomic interaction renormalizes 
the frequency of the dipole transition Cu oc c^of , whereas oc lj^ and / oc u^. 

The right hand side of Eq. ^Id^ can be conveniently rewritten as a sum of the coherent 
{i 7^ n) and incoherent {i = n) parts 

1(5^) = -7o((l + 281) tl) - 70 Y^^nSfS- + S+^rff), (21) 
where we used the identities S^S^ = ^ + Si and [S*^, f „] = 0. 
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In order to decouple correlations in Eq. (|T^. we use the Bloch representation for the 
wave functions of two- level atoms |$) js^ 

N 

K ^j) = sin |e-^-^-/^| T), + cos |e^-^-/^| j), (22) 

and obtain the closed system of equations for the averages of atomic variables {Sp = 
— |cos6'j, (S^) = |sin6'je^"^J of the following form: 

^{S'„) = - (1 + 2(5:)) (f^) - {{S:^i){S-) + (5+)(5-^i)) (Kn+i) - 

- {{sti){s-) + {s:){s-_,)){k^^,)- 



A. 



+ 2 (S'n){St){Tf), (23) 



where r = •jot, a = cjo/to ^ 1 and 

{kn±i) = 1 + /3 (3 + 3/3 + /52) Q - (5:±2)) , (25) 
(K±i) = (fLi^:) = (5:)-/9(3 + /5^)fU(5:±2)raV¥(l + (^n±2))- (26) 



There is a rapidly oscillating term ia{tn) (S^) in the system of equations (j2Sl)- In contrast 
to the case of non- interacting atoms with Fj = 1, this term cannot be removed by the phase 
shift {Sf) e^*"^(5'^). In Appendix IHl we apply the method of multitime scales Q,Q to 
eliminate the rapidly oscillating terms and deduce the following equations for slowly varying 
atomic amplitudes 



<ir) = i{S^)ir,r'))^„ a;^(r) = (5±)(r, r') exp / dt"(r„)(r,t") , (27) 
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where (...),-' denotes averaging over the time r' = ar of fast motion, 



where 



and functions 



e 



;27r(ri-rj) _ 



Wij = -I 



27r(r, - r, 



(28) 



r„ = 1 - /3 + , (29) 

Tl = l-(3{?> + (3') + + \(3' (1 + 4<^i<_i) , (30) 

= 1 + /5 (3 + 3/5 + /32) Q - ^^2) , (31) 

En±i = < - (3 + /3^) (1 + 4<±2<) + (1 + <i2) (32) 

describe correlations due to the direct interatomic interaction. 



A. Order parameter and superradiance transition 



It is known |21| that relaxation of the atomic subsystem essentially depends on the number 
of atoms A^. For non-interacting atoms ( J = and Fj = 1), the first term on the right hand 
side of Eq. fET|) describes processes of spontaneous incoherent decay. 

Transition to the regime of superradiant relaxation, where correlations are caused by 
interaction with the electromagnetic field in the vacuum state, occurs when the term de- 
scribing the coherent part in Eq. (pT|l (a sum with i ^ n) start to play the dominating role. 
Such transition can be conveniently described using the order parameter of the form (see, 
e.g. 0): 



and takes place if the number of atoms exceeds its critical value Nc 
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The presence of the factor f j = 1 — /? {S^^i + S-_i) in Eq. (jHHjl indicates that the direct 
interaction introduces additional correlations into both the coherent and the incoherent parts 
of radiation. It is the effect of such correlations on the incoherent part is responsible for the 
regime of pulse radiation js^ discussed in Sec. HI 

In Fig. 01 we present the results for dependence of the order parameter 



C^C = lim - 



dt- 



(34) 



Ef(i + 2af)r3 

computed by solving the system ()28|1 numerically at different values of the interaction param- 
eter (3 = J/hhjQ. The curves show that the critical number of atoms and the supperradiance 
transition are almost insensible to the interatomic interaction. So, similar to the case of non- 
interacting atoms, we have two qualitatively different regimes of relaxation: the incoherent 
regime at < A^^c and the regime of superradiance at > A^^c. 

0.5 

C 



0.25 




40 TV ^0 



FIG. 3: Order parameter C defined in Eq. (|34() as a function of the number of atoms N. The 
transition to the superradiant regime is shown for three different values of the ratio /3 = J/Hujq: 
(a) /3 = 0, (b) /3 = 0.5, and (c) /? = 0.9. 



B. Collective pulse relaxation 



Let us consider decay of initially inverted system of interacting atoms at subcritical 
number of the atoms A^ < A"c in more details. Dependence of the relaxation rate on time 
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calculated by solving Eq. (j^HJ with the initial conditions {S^){0) = 1/2 is shown in Fig. HJ 
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FIG. 4: Relaxation rate 7(t) defined in Eq. (|2n|) as a function of time for 10 interacting atoms at 
different values of the interaction parameter /? = J/hujQ. Solid lines represent the numerical results 
computed by solving the system (fT7|l . Dashed line is computed from the solution of Eq. (|H7|l . 

Referring to Fig. |3J it is seen that, in the absence of interaction (J = and Fj = 1), 
relaxation occurs as an incoherent spontaneous decay of excited atoms. This regime is 
characterized by a monotonic exponential decrease of the relaxation rate. 

When the interaction parameter j3 increases, spontaneous processes are suppressed and 
the time dependence of the relaxation rate reveals a non-monotonic behavior with a pro- 
nounced peak (see Fig. IH). This is the regime of collective pulse relaxation induced by the 
direct interatomic interaction. 

Collective pulse relaxation can be easily described using the mean field approximation 
for decoupling of correlations. To this end we assume that, for small number of atoms, 
the relaxation rate (j?H) is predominantly determined by the first term (incoherent part) 
describing incoherent mechanism of decay. So, the equation (1211) takes the simplified form: 

iiSn) ^ -7o((l + 2^:)f^) ^ -70 (1 + 2(^:)) (f^) (35) 
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We can now substitute the ansatz 



(^:) = (^^) + 5(^:), (36) 

assuming that {S^) ^ ^{3^), into Eq. and use the relation (IHUj) to derive the equation 
governing dynamics of atomic subsystem 

^ = -7o(l + 2{S^)){T'), {S^m = 1 (37) 

where 

(r^) ^ 1 - 2/5 (3 + P') {S') + (l + 4(^^)2) . (38) 

It can be shown that the mean field solution {S^){t) is a step-like function of time and 
is stable at \{S^) \ ^ |(5(S'^)|. The curves calculated in the mean field approximation for the 
relaxation rate 

Tl 

are shown in Fig.|3]as dashed lines. It can be seen that the results are in excellent agreement 
with the data of numerical analysis. 

It should be noted that the peak intensity is proportional to the number of atoms 
which is a consequence of a short-range character of interatomic interaction. By contrast, 
for the regime of superradiance, the intensity is typically proportional to A^^. 

For long-range interaction, the dependence of the peak intensity may change. We 
demonstrate this effect for the interaction of the form = hooo J2i ~ J^iji^tj ^i^j- 

In this case the expression (jH)) for the operator Fj has to be replaced by the relation 
Ti = 1 — P J2j jj^i that gives an additional N dependent factor. The modified expression 
for the relaxation rate is given by 

lit) ^ = 7oiV(l + 2{S^)){l + ^{N - l)/52 - (AT - l)[3/5 + (AT - 

+ 3/?2(Ar - - 2)(S")^ - p\N - \)[N - - 3)(S")3}. (39) 

Since the relaxation rate 7(t) reaches the peak at negligibly small (5*^), (5^) ~ 0, from 
Eq. pUj) its intensity can be estimated to be an increasing function of P proportional to A^^: 

3 



7max ~ 70^ 



i + r(iv-i)/32 



oc N\ 
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So, we have shown that the direct Ising-type interaction has a synchronizing effect on the 
system behavior. The resuh is that the regime of incoherent spontaneous decay changes to 
the regime of collective pulse relaxation characterized by an increase in the radiation time. 
In the subsequent section we discuss similar effects for the regime of superradiance. 

C. Enhancement of superradiance 

As it was previously discussed in Sec. IIV A| the effect of superradiance dominates the 
regime of relaxation of the atomic subsystem at supercritical values of the number of atoms, 

> Nc- Figure El presents the time dependence of the relaxation rate for A^ = 100 at 
various values of the interaction parameter j3. 

6000 

7(t) 
70 



3000 





0.5 1 70* 1.5 

FIG. 5: Relaxation rate 7(t) defined in Eq. 1)201) as a function of time for 100 interacting atoms at 
different values of the interaction parameter /? = J/huJo {2jq is the spontaneous decay rate for an 
isolated atom). The curves are computed by solving the system (jl?)) numerically. 

At /5 = 0, the peak with the intensity proportional to A^^ corresponds to the superradiance 
effect in the system of non-interacting atoms. In this system interatomic correlations are 
solely due to interaction between atoms and electromagnetic field. 

Synchronizing effect of the Ising-type interaction characterized by the parameter f3 man- 
ifests itself in enhancement of superradiance and an increase in the peak intensity. 
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The effect of enhancement can be analyzed in the mean field approximation by applying 
the procedure described in the previous section to the coherent part (the second term) of 
the relation (I2H)- So, the equation for the average of atomic population (S^) is now given 
by 

^--7oivQ-2(5^)^)(r^), (40) 

where (F^) is given in Eq. ()38|1 and the identity (the conservation law for the pseudo-spin) 
i {S+S' + S'S+) + {S'f = I is used. 

The peak of the relaxation rate 7(t) is at (S^) ~ and Eq. ()4()|1 combined with the 
relation (jnH|) provide an estimate for the gain factor (r^)!^^^^^^ oc (1 + |/3^). As for the case 
of N < Nc, the interatomic interaction delays the time of radiation. 

For long-range interaction, following the lines of reasoning presented in Sec. IIVB[ we 
derive dependence of the relaxation rate on the number of atoms N 

3 r 3 

7max ~7WI(s->=o = 7oA^^2 



1 + -{N-I)f3' 



oc A^^ (41) 



This result differs from the A^ dependence given at the end of the previous section for the 
regime of collective pulse relaxation. 

V. REGIME OF MULTIPHOTON RELAXATION 

It is known that interaction between atoms at interatomic spacin gs s maller than the wave- 
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43l . At weak interaction 



length may bring about multiphoton processes in atomic system 
with J < hujo, probability of multiphoton transitions is much smaller than that of single- 
photon transitions. So, the first order of the perturbative expansion over the atom-field 
interaction can be used to describe dynamical behavior of the system. At strong interac- 
tion with J > hwQ, single-photon processes are non- resonant and multi-photon transitions 
start to play an increasingly important role. This is the so-called regime of multiphoton 
relaxation. 

In order to illustrate how this regime may occur we qualitatively consider a model system 
of two two-level atoms with the Hamiltonian Ha = huo {Sf + Sf) — JS'fS'f. Its energy 
spectrum is schematically represented in Fig. IHl 

If the Hamiltonian of atom-field interaction has the form of a sum Hint = -f^+e*'^* + 
i^_e~*^*, where if+ {H-) describes photon emission (absorption), then the second order 
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(a): hcuo > J/2 



(b): hujo < J/2 



U) 



TT) 



IT) 



FIG. 6: Energy spectrum of two interacting atoms with Ha = + 5*1) — JSfS2- 

perturbative expression for the probabiUty of two-photon transition |0)| TT) ^ |2)| ii) is 
given by 

(ii \H^\ Ti)(Ti \H^\ TT) , (ii l^+l iT)(iT \H+\ TT) ' 



6{2uJo-2uj). (42) 



hiuoil- l3/2)-nLU huJo{l-/3/2)-fku 
For strong interaction with J > Huq, the energy level of one-particle excited states, | Ti) 
and I iT), is higher than the energy of the state | TT) (see Fig. EJo), so that one-photon 
relaxation of the two-particle excited state is excluded by the energy conservation law. As 
a consequence, relaxation occurs as a two-photon process. 

Probability of multiphoton transitions for excited atoms in free space is small because 
its magnitude is determined by A^-th order term of perturbative expansion over the atom- 
field coupling constant gks,i- So, we arrive at the conclusion that the effect of interatomic 
interaction can be the formation of long living excited states in the system. 

This conclusion, however, is not strictly valid if the density of electromagnetic modes 
has a singularity near the energy of atomic transition, hujQ. An example is a high quality 
resonant cavity. In this case, dynamical behavior of the atom-field system is characterized 
by Rabi oscillations involving many-particle and multiphoton states (44l |. 

We demonstrate this effect for the two atomic system embedded into a high-Q single-mode 
resonant cavity. The Hamiltonian of the model is 



H = hiooa+a + huoo {S'^ + S^) - JS'^S^ + kg [a (5+ + ^2+) + a+ (5f + S^)] 



(43) 



where a{a^) is the photon annihilation (creation) operator for the cavity mode, g is the 
coupling constant of the interaction between atoms and cavity mode. 
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The model ()43p is exactly solvable and the wave function for the system initially prepared 
in the state \ip{0)) = \ tt? where n is the number of photons, can be written in the explicit 
form: 



\m) 



-iuJo{n+l)t ^ 



2n + 3 
+ V(n + l)(n + 2) 



n + 2 -J,. o'^ + l 



ADt 



ADt 



D-J' 

iDt 



+ 



+ 



D + J' 

AJ't 



D \D-J' D + J' 



smDt{\ Itn + l) + I Ti,n+ 1)) 



2n + 3 



U,n + 2)- 



(44) 



where 



J' = J /Ah, D = ^/{J'y + 2g\2n + 3). 

In the absence of the interatomic interaction (J = 0), the amplitudes of one- and two- 
photon processes are of the same order of magnitude. But, at strong interaction with 
J ^ hg^/n, the two-particle two-photon amplitudes dominate and the wave function 
can be approximated as follows 

n+1 n + 2 



-i{ujo{n+l)-J'-A)t 



2n + 3 



-iAt 



2n + 3 



TT,n) + 



^. v/(n + l)(n + 2) 
+ — — smAt| ll,n + 2} 



2n + 3 



(45) 



where 



and for n ^ 1 



- (2n + 3) 



\m) 



-i{u>on— J')t 



2 J 



cosAtI tT,r^) +isinAt| ii,n + 2)). 



This means the build-up of two-atomic two-photon Rabi oscillations and can be used to 
generate a Schrodinger cat-like entangled atom-field state. 

Thus, interparticle interaction may give rise to the regime of many-particle multiphoton 
dynamics. When interaction is strong as compared with the energy of atomic transition, 
inhibition of one-photon processes in the system is accompanied by transition from the 
regime of cooperative pulse radiation (superradiance) to the generation of Fock state of light. 
In a cavity, multiphoton dynamical effects come into play under the condition J ^ hg^/n, so 
that the level of intermediate energy is essentially shifted away from one-photon resonance 
(see Fig. ini). 



23 



VI. DISCUSSION AND CONCLUSIONS 



We studied effects of direct interatomic interaction on collective processes in atom-photon 
dynamics using, as an example, a simple model of two-level atoms with Ising interaction of 
ferromagnetic type. We have found that this interaction influences radiation processes of 
atomic ensemble acting as an additional synchronizing factor. 

For weakly interacting atoms at J < hwQ, we have shown that interatomic interaction 
results in inhibition of incoherent spontaneous decay of atoms and dynamical behavior of 
the system is governed by the regime of collective pulse relaxation. This regime, though 
it bears a resemblance to superradiance, has nothing to do with the effect of phase syn- 
chronization induced by fluctuations of the electromagnetic field. For a example, in solid 
state structures, collective pulse relaxation is caused by inelastic exciton-exciton scatter- 
ing and is characterized by quadratic dependence of the radiation peak on the number of 



particles 



4G 



We have also found that interaction induced synchronization enhances superradiance 
and can be responsible for anomalous dependence of the radiation peak on the number of 
particles. In the presence of interparticle interaction collective pulse radiation and enhanced 
superradiance are both characterized by an increase in the delay time of emission. 



At the end of Sec. Ill Bl we have pointed out that when the excited state is not fully 
inverted its relaxation can be determined by transitions with frequencies equal to the atomic 
resonant frequency uq. If we consider non-excited atoms of the initial state as defects of 
the atomic chain, such relaxation scheme can be referred to as the solitonic mechanism of 
relaxation. The solitonic mechanism implies that in the course of relaxation atoms undergo 
transition to the ground state successively one after another. So, such behavior can be 
interpreted as a defect motion. More generally, since the system of equations (j^Hj) governing 
dynamics of atoms at J < hujQ is similar in structure to the Volterra system and the Toda 
lattice it might be expected that the system posses soliton-like solutions. 

Figure m presents the numerical results for the excited state with one initially non-excited 
atom. The number of atoms is small, N < Nc, and the relaxation rate is computed from 
the mean field equation (jHTj). 

Referring to Fig. [71 it is seen that relaxation of the atomic ring can be described as motion 
of the defect (Bloch wall). When the interaction parameter (3 approaches the limit of strong 
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FIG. 7: (a) Relaxation rate 7(t) as a function of time at different values of the interaction parameter 
(3. (b) Averaged population of ith atom (S"?) as a function of time and atomic number i sX (3 = 0.99. 
The curves are computed by numerically solving Eq. (|37j) for a ring of 20 atoms. One atom is 
initially in the ground state representing a point defect in the atomic ring. At sufficiently strong 
interaction, relaxation of the atomic ring clearly demonstrate the solitonic mechanism when atoms 
undergo transition to the ground state consecutively one after another. 

interaction, the relaxation rate 7(t) assumes the kink-like form and is determined by defect 
velocity. The peak is caused by collective pulse relaxation indicating that the retardation 
time of emission is too short for the defect to travel through the ring. 

However, it should be noted that we have neglected the dipole-dipole interaction by 
considering permutationally invariant atomic configurations, the approximation that is not 
valid in the presence of defect in the system. This interaction may have a destructive effect 
on the solitonic mechanism of relaxation. 

By contrast to the case of weakly interacting atoms where dynamics of the atomic sub- 
system is governed by one-particle one-photon transitions, at strong interatomic interaction 
with J ^ fiujQi these transitions appear to be excluded. In this case multiphoton transitions 
will determine relaxation of the excited atomic subsystem. This is what we called the regime 
of multiphoton relaxation. 

The regime is characterized by transition from generation of superradiant pulse to gen- 
eration of Fock quantum state of light. We have used a simple model of two atoms in a 
high-Q single mode cavity to show that such transition is accompanied by Rabi oscillations 
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involving many-atom multi pho ton states. In other words, it means generating many-particle 
entangled atom- field state 3]. 

Interestingly, transition to multiphoton dynamics is analogous to the Mott-insulator 
quantum phase transition in optically trapped atomic systems where Fock state charac- 
terizing the number of localized atoms is formed Q| . In our case, the energy of atomic 
transition Hujq and the interatomic coupling J play the role of kinetic and potential energy, 
respectively. So, the transition to generation of the Fock state of light takes place when the 
potential energy J becomes greater than hiVQ. 

In lattice atomic systems with inelastic tunneling transitions between neighboring wells, 
the Mott or Peierls transitions and the transition to multiphoton relaxation can be related 
to each other. We illustrate such a possibility by the simple example of an one-dimensional 
periodic chain of potential wells. Fig. |H1 




J 



IT), U> 



J 



■ \i>i) 



IT) 
li) 



FIG. 8: (a) Interacting particles embedded into potential lattice in the case of half-filling. Ground 
state is double degenerate corresponding to the particles localized in either odd or even wells 
(I T) = I TT ••• T) and I I) = I il ... D). Tunneling transitions lift degeneracy so that two lowest 
levels correspond to the many-particles entangled states |^o,i) ~ ^) ^\ i))- Excited states 
are separated by the gap of the width oc J. (b) For energetically non-equivalent wells, tunneling 
transitions are inelastic. For J >> fiiOQ, they can be accompanied by multiphoton processes. 



In the half-filling case, when there are half as many atoms as there are wells and the 
energy of repulsion between neighboring particles J is much larger then the hopping energy 
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of tunneling between neighboring wells Et, the ground state for identical wells should be 
degenerate in energy, which corresponds to the atoms localized in either odd or even wells, 
see Fig. |Hk- Such ordering of particles is similar to the Wigner crystal or the charge density 
waves observed in low- dimensional conductors Possibility of the density waves Mott- 

insulator phase for atoms embedded into an optical lattice was also discussed in Ref. |25|. 
These ordered states are separated from the excited ones by the gap of width ~ J. 

Local tunneling transitions may lift this degeneracy of the ground state. In this case, at 
J > Et, the ground state can be described, at least for a short chain, as a many-particle entan- 
gled state of the Schrodinger cat type, i.e., ~ (I TT ••• T) i I ii ••• i)) (the so-called GHZ 
state 6n|), associated with many-particle tunneling oscillations between odd and even wells, 

n n 

(see, e.g., Ref. which is similar to the tunnehng creep of a charge density waves |64| . 

If the wells are energetically inequivalent, Fig.|H|D, the tunneling many-particle transitions 
from the excited state |t) = lTT---T) ^ire inelastic and can be accompanied by multiphoton 
processes or the cooperative dynamic regimes discussed above. 

But our model is oversimplified and, strictly speaking, cannot be applied to inelastic 
tunneling transitions of atoms between energetically different wells. The Hamiltonian of 
atoms Ha need to be modified to take into account overlapping of the particle wave function 
in the neighboring well. This requires additional terms proportional to S^. Symmetry of 
interaction between particles and field E also has to be of more complicated form Hint oc 

Another limitation of our model is neglecting multipole transitions that, in the case of 

nn 

strong interaction, may compete with multiphoton transitions |13L I5fll|. 

Among important omissions in this paper are spectral characteristics of the radiation. In 
our case there are only two modes with the frequencies Jtuq ± J. We also have not discussed 
superradiant and subradiant Dicke states. These have the standard form because the Ising- 
type interparticle interaction makes the atomic spectrum non-equidistant but does not affect 
symmetry of these states. 
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APPENDIX A: THE BORN-MARKOV APPROXIMATION 



In this section we concentrate on the case of weakly interacting atoms with J < hujQ and 
describe the dissipative dynamics of the atomic subsystem using the Born approximation. 
To this end we apply to the master equation (jHI) the standard technique of elimination of 
bosonic variables 52j |. 

For simplicity, we consider the zero temperature case and choose an initial state with no 
correlations between the atomic subsystem and the vacuum field. So, the density operator 
can be taken in the factorized form 

Substituting Eq. (jH)) into the equation (jHl) and taking the trace over the field variables, 
we obtain 



+ [Rj{t),p{t - r)5re-^'^°r"'(*--)]) + H.c.}, (Al) 
where the sum over is transformed into the integral form 

T7" POO PTT pItT 

and the functions Rj(t), Fij{uj), p{t), r] are described by Eqs. (fTT|l - (fT3|l . 

The next step is to employ the Markov approximation. This approximation implies 
that a reservoir relaxation time is much shorter than a time-scale of the atomic subsystem 
evolution, so that p(t — r) can be replaced by p(t) in the right hand side of Eq. ()A1|) and 
the upper limit of the integral can be extended to infinity. 

For J < hujQ, the eigenvalues of the operators Fj are all positive. Using the relation 

lim / rfre^*^" = 7r5(x) ±2-, (A3) 

t^oo Jq X 

where the symbol V stands for the Cauchy principal value of the integral, we can perform 
the integral in Eq. (jAlj) and derive Eq. (fTTH) . We have used the operator functions defined 
on the eigenvalues of Fj so that iiVi\%lj) = Ti\ip) then, e.g., 6{uj — uJoTi)\ip) = 6{uj — uoTi)\ip) . 

This result can also be obtained directly by representing the terms exp(±za;of ^r) that 
enter Eq. ()A1|) as follows 
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where 



(A5) 



The principal values of the integrals in (|T0|) with i ^ j can be estimated as: 



V 



1 - (d ■ f,,f 



1 - 3(d ■ r,,; 



[l-3(d-r,,)^]i?(A:of,r,,)-[l-(d-r,,)^] 
{kQtiTijY 



A{kQt.irij) + 



(A6a) 



V 



[1 -3(d-r 



(kotirijy 



(d ■ r-.,)" 



cos{kotiri 



+ 



cos(fcofirij) 
(A;ofirij)3 



where /cq = ojq/c, 



(A6b) 



y4(x) = sin(x) ci(x) — cos(a;) si(x), 

B{x) = sin(x) si(x) + cos(a;) ci(x), x > 0, 

and si(x) = £ '^dt and ci(x) = ^dt 

In the limit VijUj/c — > 1, Eq. ()A6|1 assumes the asymptotical form p4|l coincident with the 
standard expression for the dipole-dipole interaction that does not depend on the coupling 
constant of Ising interaction J. 



APPENDIX B: TWO INTERACTING ATOMS: EXACT SOLUTION OF THE 
MASTER EQUATION 

For the case of two non-interacting atoms with J = 0, exact solution of the master 
equation (|T7j) was previously obtained in Ref. jbjj. It was shown that the decay of initially 
exited atoms demonstrates the superradiant regime and the dipole-dipole interaction does 
not influence the cooperative behavior of atoms. 

In this section we show that the master equation ()17|) for two Ising-like interacting atoms, 
J 7^ 0, can be solved along similar lines. 

The Hamiltonian of two interacting atoms is given by 

Ha = hujo {S^ + S^) - JS'^S^. 
29 
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In the basis of atomic states | H), | t|), | |t), | ||) the density operator p can be written 
the matrix form 

/pii o\ 

P22 P23 
P32 P33 
VO P44/ 

and Eq. (fTTjl reduces to the system 

Pii = -4mpii, 

Xo = 4mpii — 2n(xo + Xi), 

Xi = 4mpii — 2n(xi + a;o), (B2) 
X2 = —2nx2 + i2Qx3, 
X3 = i2Qx2 — 2nx3, 
P44 = 2n(xo + xi), 

where xq = P22 + P33, xi = P23 + P32, X2 = P23 - P32, X3 = P22 - P33, n = 70(1 + P/2f, 
m = 7o(l — /3/2)^ and Qu = ^21 = ^ is the constant of the dipole-dipole interaction. 

For the initial conditions pii(O) = 1, P22(0) = P23(0) = P32(0) = P33(0) = P44(0) = 
describing two initially excited atoms, it is easy to write down the solution of Eq. ()B2|1 

-Amt 

m 



Pii = e 

Xq = Xi 



P44 = 1 + 



n — m 
1 



-Amt 



-Ant 



) 



(B3) 



(ne-*"' - me-""' 



) 



m — n 
and also X2 = X3 = 0. 

It is seen that similar to the case of non-interacting atoms the solution ()B3|) and the 
decay rate of the atomic subsystem 

d 



lit) 



dt 



Tr{p{St + S^)}=pu-pii 



(B4) 



are both independent of the constant of the dipole-dipole interaction Q. 



APPENDIX C: ELIMINATION OF RAPIDLY OSCILLATING VARIABLES 



In order to eliminate frorn (1231) the rapidly oscillating terms m(f„)(S'^) we apply the 



method of multitime scales |56l. Is?] representing the atomic variables {S'n) as functions of 
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two time scales (5'„)(r, r'), where r' = ar {a ^ 1) is a characteristic time of fast motion. 
We use the following power series expansion of (S^) and (S^) over the small parameter 

{S^){r,r') = 5^(r')a^(r,rO = 5^(r') (^a^(r) + ^<i(r)7^(r') + , (Cla) 

{S'n){r, r') = r^T, r') = a^^^r) + ^al^[r)Y^{r') + . . . , (Clb) 

a 

where 7^(t'), 7^(t') and (5^(t') are rapidly oscillating periodic functions, so that 

Taking into account that 
and substituting the expansion (jC^laj) into Eq. (f^ . we deduce the following relation for cr+ 

+2((E„+,)a+^iC+i + (^n-i)a+_iei)+2 (^:)(f?)a+5+ (C3) 

In order to eliminate the imaginary term in the right hand side of Eq. ()C3|) we choose 
(5^(r') in the form 

C(r')=exp|z / (f„)(r,Orfd ^e^^-' 



/o 

where r„ is the zero order term of the expansion for (f„)(r, r') given by ()29|). to yield the 
equation 

Averaging Eq. over r' and retaining only the lowest order of the correlations, we 

have 

= -r^cr^ + 2 (-E„+io-.+ _iU7„+i,„ + +2 ^ r^o-+cr^Wi,„, (C5) 

where Wij is given by 

gi27r(ri-rj) _ ]^ 
Wij - 



^27r(r, - r,) ■ 
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Equation ()C5|) is the second equation of the system (PHj). The complex conjugate of 
Eq. ()C5|) gives the equation for cr^. Along the same line averaging the equation for (S^) 
from Eq. provides the last equation of Eq. (j2HI)- 
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